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Abstract 

We study dynamo action in a convective layer of electrically-conducting, compressible 
fluid, rotating about the vertical axis. At the upper and lower bounding surfaces, perfectly- 
conducting boundary conditions are adopted for the magnetic field. Two different levels of 
thermal stratification are considered. If the magnetic diffusivity is sufficiently small, the con- 
vection acts as a small-scale dynamo. Using a definition for the magnetic Reynolds number 
Rm that is based upon the horizontal integral scale and the horizontally-averaged velocity at 
C^") the mid-layer of the domain, we find that rotation tends to reduce the critical value of Rm 

above which dynamo action is observed. Increasing the level of thermal stratification within 
the layer does not significantly alter the critical value of Rm in the rotating calculations, but 
it does lead to a reduction in this critical value in the non-rotating cases. At the highest 
computationally-accessible values of the magnetic Reynolds number, the saturation levels of 
the dynamo are similar in all cases, with the mean magnetic energy density somewhere be- 
O I tween 4 and 9% of the mean kinetic energy density. To gain further insights into the differences 

between rotating and non-rotating convection, we quantify the stretching properties of each 
flow by measuring Lyapunov exponents. Away from the boundaries, the rate of stretching due 
to the flow is much less dependent upon depth in the rotating cases than it is in the corre- 
ct sponding non-rotating calculations. It is also shown that the effects of rotation significantly 
1 1 reduce the magnetic energy dissipation in the lower part of the layer. We also investigate 
certain aspects of the saturation mechanism of the dynamo. 

> 

^ 1 Introduction 

m 

In a hydromagnetic dynamo, motions in an electrically-conducting fluid lead to the ampli- 
fication of a seed magnetic field. Dynamo action can only occur if the inductive effects of 
the fluid motions outweigh the dissipative effects of magnetic diffusion. Provided that the 
I magnetic energy density of the seed field is very much less than the kinetic energy density of 

the flow, the early stages of evolution of the dynamo process are effectively kinematic, which 
implies that the magnetic energy in the system grows exponentially with time (although the 
m magnetic energy will tend to fluctuate about this exponential trajectory if the velocity field 

is time-dependent). Eventually, however, the nonlinear feedback of the magnetic field upon 
the flow (via the Lorentz force) becomes dynamically significant. This halts the growth of the 
dynamo-generated magnetic field, leading to a saturated nonlinear dynamo. There are many 
examples of natural dynamos. For example, it is generally accepted that the 11 year solar mag- 



o 



netic cycle (see, for example, Stix 2004), is driven by an oscillatory dynamo. Similar dynamo 
mechanisms drive cyclic magnetic activity on other stars ( |Baliunas et al.\ |1996[ ). Dynamo 



action is also responsible for sustaining the Earth's magnetic field (see, for example, Roberts 



& Glatzmaier 2000[ |. It is also possible to drive hydromagnetic dynamos in laboratories, as 



demonstrated by recent liquid metal experiments (such as the French VKS experiment, see 



for example, Monchaux et al. 20091 



Convection plays a crucial role in many natural dynamos, particularly in the context of 
solar, stellar and planetary dynamos. Although it is possible to study global dynamo mod- 



els numerically (Brun et al. 2004|, most of the theoretical work on this problem has been 
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based upon local models of convection, where a fluid layer is heated from below and cooled 
from above. In the context of Boussinesq convection, this classical setup has been studied 
numerically by Meneguzzi & Pouquet (J1989) and Cattaneo (19991. In highly-conducting flu- 



ids, Boussinesq convection can act as an efficient dynamo, producing small-scale, intermittent 
magnetic fields. In the Boussinesq approximation, the effects of compressibility are neglected. 
More recent studies have focused upon dynamo action in local models of convection in fully 
compressible fluids ( jVogler fc Schiissler| [^507} |Brummell etaT] [2010} |Bushby et al\ [20lo] ) . It 
is now well-established that these compressible models can also drive small-scale dynamos. 

Previous calculations have clearly demonstrated that rotation is not a necessary ingredient 
for small-scale dynamo action in convectively-driven flows. However, rotation is present in most 
natural dynamos, and this additional physical effect can profoundly influence the behaviour 
of a convective layer. For example, when the rotation vector is aligned with the vertical axis, 
rapid rotation not only inhibits convection, but it also reduces the preferred horizontal scale 

Even far from onset, 



of the convective instability ( Chandrasekhar 
convection is 



1961). 



when compressible 

fully turbulent, rotation can also play an important role: helical convective 
plumes tend to become aligned with the rotation vector. This is particularly apparent when 
the axis of rotation is tilted away from the vertical ( Brummell et al. 1996 199861. Given these 



hydrodynamical considerations, we might reasonably expect the dynamo properties of rotating 
convection to differ from the equivalent non-rotating case. There have been numerous studies 



of this problem in Boussinesq convection (Childress & Soward 1972 St. Pierre 1993 Jones 



& Roberts 2000 Stellmach & Hansen 2004 Cattaneo & Hughes 2006 1 , but the compressible 



case is less well understood. Existing studies have generally adopted complex models, with 
multiple polytropic layers, inclined rotation vectors, or additional physical effects such as an 



imposed velocity shear (see, for example, Brandenburg et al. 1996 Kapyla, et al. 2008 20091. 
In fact, the simpler problem of dynamo action in a single layer of turbulent compressible 
convection, rotating about the vertical axis, is still not fully understood. Therefore, the aim 
of this paper is to study the ways in which rotation and compressibility influence the dynamo 
properties of convection in a simple polytropic layer. 

The paper is organised as follows. The governing equations, boundary conditions and 
parameters, together with the numerical methods are discussed in the next section. Consid- 
erations about the dimensionless numbers of interest, and the Lyapunov exponents of hydro- 
dynamic convection are presented in Section [3] In Section |4j we discuss results from a set of 
dynamo calculations. Finally, in Section |5j we present our conclusions. 



2 Model and method 

2.1 Model and governing equations 

We consider the evolution of a plane-parallel layer of compressible fluid, bounded above and 
below by two impenetrable, stress-free walls, a distance d apart. The upper and lower bound- 
aries are held at fixed temperatures, To and To + AT respectively. Taking AT > implies 
that this layer is heated from below. The geometry of this layer is defined by a Cartesian 
grid, with x and y corresponding to the horizontal coordinates. The z-axis points vertically 
downwards, parallel to the constant gravitational acceleration g = ge z . The layer is rotating 
about the z-axis, with a constant angular velocity ft — ttz. The horizontal size of the fluid 
domain is defined by the aspect ratio A so that the fluid occupies the domain < z < d 
and < x,y < Ad. The physical properties of the fluid, namely the specific heats c p and c„, 
the dynamic viscosity p, the thermal conductivity K, the magnetic permeability po and the 
magnetic diffusivity 77, are assumed to be constant. The model is identical to that used by 



Matthews et al. (19951, except for the addition of rotation. 



It is convenient to introduce dimensionless variables, so we adopt the scalings described 



in Bushby et al. (20081. Lengths are scaled with the depth of the layer d. The temperature 
T and the density p are scaled with their values at the upper surface, To and po respectively. 
The velocity u is scaled with the isothermal sound speed \/R,Tq at the top of the layer, where 
R, = c p — Cv is the gas constant. We adopt the same scaling for the Alfven speed, which 
implies that the magnetic field B is scaled with y/popoR*To- Finally, we scale time by an 
acoustic time scale d/V^?*T>. 

We now express the governing equations in terms of these dimensionless variables. The 
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equation for conservation of mass is given by 

| = -V-(pn). (1) 
Similarly, the dimensionless momentum equation can be written in the following form, 

p + K,oTa /2 z xu) =-v(p+ + 9 (m + 1) pz + V ■ (BB - puu + k<jt) , (2) 



dt u J V 2 . 

where P is the pressure, given by the equation of state P = pT, and t is the stress tensor 
defined by 

— ® Ui i _ 2 x du k 

Tl3 ~ <).,-, + <)., . r ij dx k ■ [6) 

Several non-dimensional parameters appear in Equation The parameter 6 = AT '/To is the 
dimensionless temperature difference across the layer, whilst m = gd/R t AT— 1 corresponds to 
the polytropic index. The dimensionless thermal diffusivity is given by k = K/dpoc p (R,To) 1 ^ 2 
and a — pc p /K is the Prandtl number. Finally, Tao is the standard Taylor number, Too = 
4p 2 ,Q, 2 d 4 / p 2 , evaluated at the upper boundary. The induction equation for the magnetic field 
is 

^=Vx(uxB - CokV x B) , (4) 

where £o = V c pPo/K is the ratio of magnetic to thermal diffusivity at the top of the layer. 
The magnetic field is solenoidal so that 

V • B = . (5) 

Finally, the heat equation is 

= -U ■ VT - ( 7 - 1)TV ■ u + ^V 2 T + ^ ~ 1} (ar 2 /2 + Co |V x B\ 2 ) , (6) 
where 7 = c p /c v . 

To complete the specification of the model, some boundary conditions must be imposed. 
In the horizontal directions, all variables are assumed to be periodic. As has already been 
described, the upper and lower boundaries are assumed to be impermeable and stress-free, 
which implies that u XiZ = u VtZ — u z = at z = (the upper boundary) and z = 1 (the 
lower boundary). Having non-dimensionalised the system, the thermal boundary conditions 
at these surfaces correspond to fixing T = 1 at z — and T = 1 + 6 at z = 1. For the magnetic 
field boundary conditions, two choices have typically been made in previous studies. One can 
consider the upper and lower boundaries to be perfect electrical conductors or one can adopt a 
vertical field boundary condition. We choose appropriate conditions for perfectly-conducting 
boundaries, which implies that B z = B XtZ — B VtZ — at z = and 2 = 1. This is partly 
to facilitate comparison with the Boussinesq study of |Cattaneo fc H ughes (2006), but this is 
not the only motivation for adopting boundary conditions that enforce B z = at the upper 
surface. Strong concentrations of vertical magnetic flux at the upper surface tend to become 
partially evacuated (see, for example, |Bushby et aL[ [2008 ) , which dramatically increases the 
local Alfven speed (as well as significantly reducing the local thermal diffusion time scale). 
This can impose very strong constraints upon the time-step in any explicit numerical scheme. 
Hence, there are also numerical reasons for adopting perfectly-conducting boundary conditions. 



In this context, it is worth noting that it has been shown in the Boussinesq case (Thelen & 



Cattaneo 2000) that the dynamo efficiency of convection is largely insensitive to the detailed 



choice of boundary conditions for the magnetic field (see also the compressible calculations 



of Kapyla, et al. 20101. Although not reported here, we have carried out a small number of 
simulations with vertical field boundary conditions, and the results are qualitatively similar 
in that case. 

2.2 Model parameters and initial conditions 

With the given boundary conditions, it is straightforward to show that these governing equa- 
tions have a simple equilibrium solution, corresponding to a hydrostatic, polytropic layer: 

T = 1 + Oz , p=(l + 6z) m , u = , B = . (7) 
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Run 


Ra 


Ta 





K 


Re 


Rl 


3 x 10 5 





3 


0.0055 


153 


R2 


4.6 x 10 5 


10 5 


3 


0.0044 


157 


R3 


5 x 10 5 





10 


0.0242 


149 


RA 


6 x 10 5 


10 5 


10 


0.02 


152 



Table 1: The set of parameters for the four different cases. Note that the (mid-layer) Taylor number 
is defined by Ta = 4(1 + / 2) 2m p^Vl 2 d A / p? . See the main part of the text for the definitions of the 
global Reynolds number, Re, and the Rayleigh number, Ra. 



This polytropic layer (coupled with a small thermal perturbation) is an appropriate initial 
condition for simulations of hydrodynamic convection. All the hydrodynamic flows that are 
considered in this paper are obtained by evolving the governing equations from this basic 
initial condition. Once a statistically-steady state has been reached, the dynamo properties of 
these flows can be investigated by inserting a weak (seed) magnetic field into the domain. 

There are many non-dimensional parameters in this system. Since it is not viable to 
complete a systematic survey of the whole of parameter space, we only vary a subset of the 
available parameters. Throughout this paper, the polytropic index is fixed at m = 1, whilst 
the ratio of specific heats is given by 7 = 5/3 (the appropriate value for a monatomic gas). 
These choices ensure that the initial polytropic layer is unstable to convective perturbations. 
We also fix the Prandtl number to be a = 1. In order to study the effects of varying the 
stratification, we adopt two different values of 9. The case of 9 = 3 corresponds to a weakly- 
stratified layer, whilst 9 = 10 represents a highly-stratified case in which the temperature and 
density both vary by an order of magnitude across the layer. 

The main aim of this study is to address the effects of rotation and compressibility upon 
dynamo action in a convective layer. So, for each value of 9, we consider a rotating and a 
non-rotating case (which implies four different cases overall). Note that some care is needed 
when comparing simulations of rotating convection at different values of 9. The Taylor number 
that appears in the governing equations, Tao, corresponds to the Taylor number at the top 
of the domain. Given the differing levels of stratification, it makes more sense to specify the 
same mid-layer Taylor number for each of the rotating cases. Since the Taylor number is 
proportional to p 2 , the mid-layer Taylor number (in the unperturbed polytrope) is defined by 
Ta = Ta (l + 9/2) 2m . Two different values of Ta are considered, Ta = (the non-rotating 
cases) and Ta = 10 5 . Further discussion regarding the depth-dependence of the Taylor number 
is given in the next section. 

Another parameter that needs to be specified is the dimensionless thermal diffusivity, re. 
However, rather than prescribing values for k, we define the mid-layer Rayleigh number 

Ra = (m + 1 - my) (1 + 9/2) 2m ~ 1 (m ± ^ , (8) 

which is inversely proportional to k 2 and measures the destabilising effect of the superadiabatic 
temperature gradient relative to the stabilising effect of (viscous and thermal) diffusive pro- 
cesses. As described in the Introduction, rotation tends to stabilise convection, whilst larger 



values of 9 also have a stabilising effect (Gough et al. 1976). Therefore it does not make sense 



to keep Ra constant in all cases. Instead, we vary the Rayleigh number from one calculation 
to the other, ranging from Ra — 3 x 10 5 up to Ra — 6 x 10 5 . The aim was to reach similar 
values of the Reynolds number in all cases. We shall discuss another possible definition for 
the Reynolds number in the next section, but for now we define a global Reynolds number, 
based upon the rms velocity (J7 rma ), the kinematic viscosity at the mid-layer (ner/pmid, where 
Pmid is the mean density at the mid-layer of the domain) and the depth of the layer (which 
equals unity in these dimensionless units). Hence this global Reynolds number is given by 

Re = pmidj7rms . (9) 

ft<7 

The choices for Ra that we have used imply that Re is approximately 150 for each of the four 
cases. A summary of our choice of parameters for each case is given in Table [l] 

All the parameters that have been discussed so far relate to the hydrodynamic properties 
of the convection. If all other non-dimensional parameters in the system are fixed, the early 
evolution of any weak magnetic field in the domain depends solely upon value of £0, which is 
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a parameter that can be set as the field is introduced. This parameter is proportional to the 
magnetic diffusivity, 77, so we require low values of Co for dynamo action. Equivalently, we could 
say that dynamo action is only expected in the high magnetic Reynolds number regime. As 
for the Reynolds number, we shall discuss an alternative definition for the magnetic Reynolds 
number in Section [3. 1| However, for the moment, we make an analogous global definition: 

Rm = ^P . (10) 
«Co 

We choose a range of values for Co which give values of Rm that vary between approximately 
30 and 480 (0.12 < Co < 2.0 for 9 = 3 and 0.05 < Co < 0.8 for 6 = 10). This range of 
values is restricted by numerical constraints: smaller values of Co (higher magnetic Reynolds 
numbers) would require a higher level of numerical resolution, which would greatly increase 
the computational expense. 



2.3 Numerical method 

Solving the equations of compressible convection in the presence of a magnetic field is more 
demanding numerically than equivalent Boussinesq calculations, so it is important to use a 
well optimised code. The given set of equations is solved using a modified version of the mixed 



pseudo-spectral/finite difference code originally described by Matthews et al. (19951. Due to 
periodicity in the horizontal direction, horizontal derivatives are computed in Fourier space 
using fast Fourier transforms. In the vertical direction, a fourth-order finite differences scheme 
is used, using an upwind stencil for the advective terms. The time-stepping is performed by an 
explicit third-order Adams-Bashforth technique, with a variable time-step. For all calculations 
presented here, the aspect ratio is A = 4. The resolution is 256 grid-points in each horizontal 
directions and 120 grid-points in the vertical direction. A poloidal-toroidal decomposition is 
used for the magnetic field in order to ensure that the field remains solenoidal. As explained 
later, we also aim to calculate Lagrangian statistics. Trajectories of fluid particles are therefore 
computed using the following equation: 

^f=u( Xp ,t), (11) 

where x p is the position of the particle. The velocity at the particle position is interpolated 
from the grid values using a sixth-order Lagrangian interpolation scheme. The boundaries 
are treated with a decentred scheme and we carefully check that all the particles remain in 



the fluid domain. The time-stepping used to solve Equation (111 is the same as for the other 
evolution equations in the system. 



3 Hydrodynamical considerations 

As we have already described, fully-developed hydrodynamic convection is used as a starting 
point for all dynamo calculations. In this section, we consider the properties of the hydrody- 
namical flows that are obtained in each of four cases that are listed in Table [T] 



3.1 Dimensionless numbers 

In a stratified layer, most quantities of interest will be a function of depth. This is true not 
only for the thermodynamic quantities in the flow, but also for some of the parameters in 
the system. In the previous section, we defined the Taylor number and the Reynolds number 
using mid-layer values for the density, in addition to using the layer depth for the characteristic 
length scale. These are certainly valid definitions for these quantities, but further insight can be 
gained by considering the depth-dependence of these parameters, particularly when comparing 
calculations with different levels of stratification. This is a point that we consider in detail in 
this subsection. 

For each value of z, it is possible to define local dimensionless numbers by considering 
horizontally-averaged quantities. Under such circumstances, it is more sensible to define these 
parameters in terms of a horizontal length-scale rather than using the depth of the layer (which 
equals unity in this dimensionless system). We choose here the horizontal integral length scale, 
defined by 

lo{z)= (u^Ux^J , (M*,t)ui( X + re x ,t)) z dr (12) 
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Integral length scale /q(z) Root mean square velocity t/rms(z) 



Figure 1: Variation of (a) the integral length scale lo(z) and (b) the root mean square velocity 
Urms(z) with depth for the four calculations considered here. The rotating cases correspond to the 
empty blue symbols whereas the non-rotating cases correspond to the red crosses. 



(where we have assumed the the flow is horizontally-isotropic). Here, and in the following, 
the brackets < . > z mean a statistical average over horizontal coordinates and time at a given 
depth z, whereas the brackets < . > (with no subscript) mean a statistical average over time 
and all spatial coordinates. We note that this horizontal scale (unlike the vertical dimension 
of the domain) will vary not only with depth, but also from one flow to another. 

Figure [I] shows the variation with depth of both the integral scale, Iq(z), and the local 
rms velocity, U T m B (z) =< \u(x,t)\ 2 >V 2 , for the four cases. Each of these horizontally- 
averaged quantities is averaged over more than 100 acoustic crossing times. In all four cases, 
the trends are similar. Near the surface, the integral scale lo(z) decreases with depth as the 
flow becomes more turbulent, and therefore less correlated. As the flow reaches the lower 
part of the layer, the integral length scale increases again, presumably due to the presence of 
boundaries. Comparing the two non-rotating cases, one observes that the effect of increasing 
9 is to decrease the integral scale of the flow. However, whatever the value of 9 is, the length 
scales are very similar in the two rotating cases, and are always smaller than the corresponding 
integral scales in the non-rotating calculations. This mirrors the result from linear theory that 



rotation tends to decrease the size of the convective cells ( Chandrasekhar 19611. The rms 
velocity is comparatively independent of depth, except close to the lower boundary where (like 
the integral scale) it increases. Note that U TUls (z) is larger in both the highly-stratified cases. 
At fixed stratification, the rms velocity is smaller in the rotating cases. 

We can now use these velocity and length scales to define local dimensionless parameters 
for this system. Unlike their corresponding global values (based upon the depth of the layer 
and the depth-averaged rms velocity), these are always functions of z. We now define the local 
Reynolds number to be 

= (P0M)),tW*)*o(*) (13) 
Similarly, we can define the local Rossby number to be 

*•« = W) ' (14) 

whilst 

4(p(x,t)) 2 z n 2 it(z) 

Ta(z) = — (15) 

gives a local definition for the Taylor number. Although we are focusing upon the hydrody- 
namic properties of convection in this section, this is also a convenient place to define the local 
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Figure 2: Dimensionless numbers versus depth. From left to right: Reynolds number Re(z), 
magnetic Reynolds number Rm{z) (the global value of Rm is 30 in each of the four cases shown 
here), Rossby number Ro(z) and Taylor number Ta{z). 



magnetic Reynolds number, 

M , )= M^M. (16) 

In the dynamo calculations in the next section, we shall often refer to the mid-layer value of 
the local magnetic Reynolds number, i.e. Rm (0.5). 

The depth-dependence of the local Reynolds number, the local magnetic Reynolds number 
(for four cases corresponding to a global value of Rm — 30), the local Rossby number and 
the local Taylor number is shown in Figure [2] First of all, it is clear that the local Reynolds 
number is an increasing function of depth, which is a clear indication that the flow is becoming 
more turbulent as z increases. This variation with depth is partly due to the fact that the 
dynamic viscosity /i is held constant in all our simulations. An important observation is that 
the local Reynolds number (and the local magnetic Reynolds number) is always smaller in 
the rotating cases. This is despite the fact that the global Reynolds numbers and magnetic 
Reynolds numbers (as defined in the previous subsection) are the same in all cases. This 
reduction in the local Reynolds number is due to the fact that rotation leads to motions with 
a smaller characteristic horizontal length scale (as well as a lower rms velocity). So, if these 
local definitions for the Reynolds numbers give a clearer indication of the level of turbulence 
in the flow, the rotating calculations are actually slightly less turbulent than suggested by the 
values of the global Reynolds numbers. We shall return to this point in Section [4] Returning 
to Figure [2] we see that the Rossby number is roughly constant across the layer, but is smaller 
in the high 9 case. Although we could conclude that this implies that the effects of the Coriolis 
force with respect to inertia are more important in the highly stratified case, the key point is 
that Ro{z) is always much less that unity, which implies that the flows are both rotationally 
dominated. Finally, the local Taylor number is largely independent of depth (and 9) in the 
middle of the domain, but increases near the boundaries. Note that this behaviour is very 
different from the global Taylor number, which actually increases with depth. This apparent 
discrepancy is due to the fact that the local Taylor number is strongly influenced by variations 
in the integral scale (varying as lg(z)), which tend to outweigh any local increases in the mean 
density. The Taylor number (in particular) highlights some of the difficulties that must be 
faced when defining appropriate dimensionless numbers in highly-stratified fluids. 

3.2 Lyapunov exponents 

From the point of view of dynamo action, it is interesting to quantify the amount of stretching 
in the flow. This can be achieved by measuring the corresponding Lyapunov exponents. 
Several previous studies have explored the relationship between dynamo action and Lyapunov 
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Figure 3: (a) Maximum Lyapunov exponents, A e , versus depth, (b) Maximum Lyapunov exponents 
scaled by the mean turn-over time lo/U Tms . Right column: Comparison between non-rotating (two 
first rows) and rotating convection (two last rows). From left to right: Temperature fluctuations, 
enstrophy and short time Lyapunov exponents (STLE) at a given depth z. From top to bottom: 
z = 0.2 without rotation, z — 0.8 without rotation, z — 0.2 with rotation, z = 0.8 with rotation. 
9 = 3 in all cases. 



exponents (see, for example, Finn & Ott 1988 Brandenburg et al. 1995 Brummell et al. 



1998 a Tanner & Hughes 20031, and we would expect regions with large Lyapunov exponents 



- log 

t 



(17) 



to be playing a dominant role in the dynamo process. Lo estimate the short-time Lyapunov 
exponent, we release 5 x 10 4 fluid particles randomly within the convective layer. For each 
particle, we also release a second particle at a distance do = 10~ 4 apart from the first particle, 
in a random direction. This pair of particles is then followed as it moves through the fluid. 
The short-time Lyapunov exponent A e (Eckhardt & Yao |1993[ ) is then calculated using the 
following expression: 

m 

do 

where d(t) is the distance between the two particles at time t. The Lyapunov exponent is 
computed after a fixed time (approximately equal to one crossing time) and the distance 
between the particles is then renormalised to do. Lhe results appear to be insensitive to 
the choice of the initial separation distance do- Furthermore, they do not appear to depend 
upon the details of the method that is used to reinitialise the particle positions. All initial 
separations will quickly align with the expanding manifold, so we would expect to see very 
little dependence upon the initial condition for sufficiently large time, t. 

To study the depth dependence of the stretching within the flow, we consider the horizontal 
average of the largest Lyapunov exponent at each value of z. In addition, we time-average the 
resulting Lyapunov exponent over approximately 20 convective turnover times. Figure |3^a) 
shows the mean maximum Lyapunov exponent versus depth for each of the four cases. In all 
cases, the stretching increases with depth (as the flow becomes more and more turbulent), 
with the Lyapunov exponents taking their maximum values at the bottom of the convective 
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domain. Unsurprisingly, this is reminiscent of the depth-dependence of the local Reynolds 
number. This depth-dependence has important implications for the dynamo problem: we 
would expect the magnetic field to be mostly generated in the lower part of the layer, where 
the rate of stretching is maximal. Compared to the equivalent non-rotating cases, in rotating 
convection we see higher values of the Lyapunov exponent at the top of the layer and lower 
values at the bottom of the layer. In other words, the maximum Lyapunov exponent is 
less depth-dependent in the rotating cases than it is in the non-rotating calculations. From 
a dynamo perspective, this suggests that we should expect the magnetic fields to be more 
homogeneously distributed across the layer in the rotating cases than they are in the non- 
rotating simulations. Figure |3|b) shows a plot in which the Lyapunov exponents are scaled 
by the local convective turnover time lo(z)/U T ms(z). We use a similar scaling for the dynamo 
growth rates in the next section, where the implications of this plot will be discussed in more 
detail. Here we simply note two key features of this scaled plot. Firstly, this scaling causes the 
four curves to collapse onto a single curve in the upper part of the layer. Secondly, we note 
that the scaled Lyapunov exponents suggest that there is less stretching in the lower part of 
the domain in the rotating cases. 

On the right-hand side of Figure [3] we show a snapshot of the temperature, enstrophy 
and the short time Lyapunov exponent (STLE) in a horizontal plane at two different depths 
[z — 0.2 and z = 0.8) for both non-rotating (top) and rotating (bottom) convection in the 
weakly-stratified (9 = 3) case. In the temperature plot, brighter contours correspond to 
warmer regions of fluid. As expected from our considerations of the horizontal integral scale, 
it is clear that the horizontal scale of convection is smaller in the rotating case than it is in 
non-rotating convection. This reduction in horizontal scale is also apparent in the enstrophy 
plot where bright regions, corresponding to areas of high (squared) vorticity, outline the edges 
of the convective cells. The STLE map is obtained by releasing 10 particle pairs at a given 
depth. The particles are followed for approximately one crossing time. The STLE is then 
computed using Equation (171 and plotted at the initial particle pair position. Comparing 
the STLE map with the enstrophy distribution, it is clear that there is a correlation between 
zones of strong stretching and regions of high vorticity. Given this correlation, it is natural 
to conclude that the higher (unsealed) Lyapunov exponent in the upper part of the rotating 
simulations, as observed in Figure [3ja), is mostly due to the larger number of convective cells 
that are present. The scaled plot that is shown in Figure^b) tends to support this conclusion: 
the scaling that has been used here takes into account this "filling factor" effect, which explains 
why all the scaled Lyapunov exponent curves collapse onto a single curve near the surface in 
that case. 



4 Dynamo simulations 

In this section, we investigate the dynamo properties of the four convective flows that are being 
considered in this paper. Each dynamo calculation is initialised by introducing a seed magnetic 
field into statistically-steady hydrodynamic convection. The initial spatial structure of this 
magnetic field is given by B(x, y, z) ={A cos(fciy), A cos(kix), 0), where fc; = 2n/X and A is the 
peak amplitude of the initial perturbation. This is almost the simplest possible magnetic field 
configuration that is consistent with the imposed boundary conditions, whilst also ensuring 
that there is no net (horizontal or vertical) magnetic flux across the computational domain. 
Test calculations suggest that the evolution of the dynamo is largely insensitive to the precise 
choice of initial conditions. 



4.1 The kinematic regime 

Initially, we consider the kinematic dynamo regime in which the seed magnetic field is assumed 
to be weak. This implies that the magnetic field tends to grow (or decay) exponentially at 
a rate that is determined by the value of magnetic Reynolds number. Given that the actual 
growth of the magnetic energy fluctuates in time, long time-averages are often needed in order 
to accurately measure growth rates. This kinematic phase of evolution can be indefinitely 
prolonged by switching off the Lorentz force terms in the momentum equation and the ohmic 
heating terms in the temperature equation. This is the procedure that is adopted in this 
subsection. For each velocity field, we carry out 6 kinematic calculations at different values 
of £()• This parameter is chosen so that the global magnetic Reynolds number ranges from 
approximately 30 to 480 in each case. This implies that the magnetic Prandtl number (the 
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ratio of the global magnetic Reynolds number to the global Reynolds number) varies from 
approximately 0.2 to 3.2. 

The main aim of these kinematic calculations is to determine the ways in which rotation 
and stratification influence the rate of growth (or decay) of the seed magnetic field. The upper 
part of Figure [4] shows two plots (for 9 = 3 and 9 = 10) of the kinematic growth rate of the 
magnetic energy, A, versus the global magnetic Reynolds number. In all cases, A has been 
scaled by the mid-layer turnover time of the turbulence, Zo(0.5)/[/ rms (0.5). This rescaling has 
been carried out to facilitate comparison between the four cases, as well as the Lyapunov 
exponent calculations from the previous section. Results from this kinematic study highlight 
several key points. Firstly, it is clear that the critical value (for the onset of dynamo action) 
of the global magnetic Reynolds number is rather similar in all four cases. It is only in the 
non-rotating 9 = 3 case that there is a suggestion of a slightly higher critical global magnetic 
Reynolds number (although the error bars in the growth rate close to criticality are large 
enough to suggest that this difference may not be significant). Adopting this definition for the 
magnetic Reynolds number, we see that the scaled growth rates in the highly stratified case 
seem to be rather insensitive to the presence of rotation. It is only in the weakly stratified 
case, at small values of the global magnetic Reynolds number, that significant differences are 
seen between the rotating and the non-rotating cases. So it is tempting to conclude from these 
results that the dynamo properties of compressible convection (particularly at higher values 
of the global magnetic Reynolds number) are largely insensitive to the presence of rotation, 
as well as variations in the level of stratification. 

However, some caution is needed when interpreting these results. As we discussed in the 
previous section, a global definition for the magnetic Reynolds number takes no account of the 
horizontal scale of convection (a quantity that varies significantly between the four cases that 
are being considered). So we would argue that it is more sensible to consider representative 
values of the local magnetic Reynolds number for the purposes of this comparison. The lower 
part of Figure [4] also shows the growth rate curves for the four different cases, but this time A 
has been plotted against the mid-layer value of the local magnetic Reynolds number, _Rm(0.5). 
It is now possible to discern clear differences between the four sets of calculations. We first 
discuss the effect of varying 9 without rotation. For 9 = 3, the critical value for Rm(0.5) is 
about 420, whereas for 9 = 10, the critical value for 7?m(0.5) is about 290. Therefore, an 
increase in the level of stratification reduces the critical value of the local magnetic Reynolds 
number. Let us now consider the two rotating cases. A striking feature of the plots in the 
lower part of Figure [4] is how similar the growth rate curves of the two rotating cases are. 
For both values of 9, the critical value of _Rm(0.5) is about 220. Therefore, independently 
of the level of stratification, rotation tends to reduce the critical value for the local magnetic 
Reynolds number. 

Whatever definition is adopted for the magnetic Reynolds number, it is clear that the 
growth rate curves exhibit certain characteristic features. The dotted lines on Figure|4j which 
are shown on these plots for indication, correspond to a logarithmic scaling (see, for example, 
Rogachevskii & Kleeorin 19971 and a square-root scaling ( Schekochihin et al. 2004| of the 



growth rate with the magnetic Reynolds number. A logarithmic scaling seems to be valid 
(at least for all magnetic Reynolds numbers considered here) for the non-rotating 9 = 3 case. 

1/2 

This situation is less clear in the highly-stratified non-rotating case, although an R AI scaling 
may be more appropriate here. The other possibility is that there are two different regimes in 
this case, with a logarithmic scaling holding only for larger values of the magnetic Reynolds 
number. The scalings are again not completely convincing in the rotating cases, but the data 

1/2 

seem to be compatible with an R^j scaling, irrespective of the value of 9. In both plots in 
the lower part of Figure [4] it is clear that the rotating cases always have higher growth rates 
than the equivalent non-rotating calculations. The largest difference in growth rates between 
the rotating and non-rotating calculations can be seen at low magnetic Reynolds number in 
the 9 = 3 case. 

The observed variations in the growth rate curves suggest that rotation is beneficial for 
dynamo action. However, as discussed in the previous section, it could be argued that the 
rotating calculations are (in some sense) less turbulent that their non-rotating counterparts: 
although the global Reynolds numbers are similar in all cases, the local Reynolds number is 
significantly smaller in the rotating calculations. So we have to consider the possibility that 
it is actually the difference in the local Reynolds numbers that is giving the impression of 
enhanced dynamo action in the rotating cases. In order to address this issue, we performed 
an additional set of kinematic simulations for a new convective flow (referred to here as case 
Rib). Setting Ta = 0, we choose 9 = 3 since this it is in this case that we observe the 
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Figure 4: Growth rates of the magnetic energy versus the magnetic Reynolds number for 6 = 3 
(left) and — 10 (right). In the upper plots, the global magnetic Reynolds number has been used. 
In the lower plots, the same growth rates have been plotted against the mid-layer value of the local 
magnetic Reynolds number, Rm (0.5). The rotating results are plotted with triangles whereas the 
non-rotating results correspond to the circles. The growth rates are normalised by the mid-layer 
turnover-time Z (0.5)/?7 rms (0.5). The error bars correspond to the maximum and minimum value 
for the growth rate when considering different time intervals. The logarithmic and square-root 
scalings are plotted for reference. 
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Figure 5: Results from an additional set of non-rotating simulations at lower local Reynolds num- 
ber, (a) The local Reynolds numbers versus depth for simulations Rl, Rib and R2. (b) Growth 
rates of the magnetic energy versus the mid-layer value of the local magnetic Reynolds number, 
Rm (0.5), for all the 8 = 3 simulations. The rotating results, from simulation R2, are plotted with 
triangles whereas the non-rotating results correspond to the circles (filled for Rl, and empty for 
R2). In all cases, the growth rates are normalised by the mid-layer turnover-time Zq(0.5) /U ims (0.5). 



greatest difference between the dynamo properties of rotating and non-rotating convection. 
All parameters are identical to case Rl (as given in Table [I]), except the Rayleigh number, 
which we reduce to Ra = 5 x 10 4 . The global Reynolds number is now significantly smaller 
than in all other calculations reported in this paper. However, the local value now has a 
similar depth dependence to the equivalent rotating case, R2, as illustrated in Figure [5ja). 
Figure [5jb) shows a comparison of the kinematic dynamo growth rates for all the 9 — 3 
calculations (Rl, Rib and R2). At similar local magnetic Reynolds numbers, the growth rates 
for the lower Rayleigh number case are slightly higher than the equivalent non-rotating higher 
Rayleigh number case (721) . In other words, the reduction in the local Reynolds number 
has increased the efficiency of the dynamo. However, there is still a significant discrepancy 
between these growth rates and those found in the rotating calculations. This indicates that 
the reduction in the critical value of J£a/(0.5) in the rotating cases cannot be explained simply 
by the differences between the local Reynolds numbers. Therefore this change in the critical 
local magnetic Reynolds number must be due to the effects of rotation. 

Before concluding this section on kinematic growth rates, it is worth commenting on the 
magnitudes of the growth rates in the dynamo regime. Obviously the range of magnetic 
Reynolds numbers is restricted by numerical constraints. However, it is worth noting that our 
positive growth rates (when properly normalised) are comparable to the values reported by 



Kapyla et al. (20081, who also considered simulations of compressible convection (albeit with 



overshoot and shear), and those found in the forced homogeneous turbulence calculations of 



Haugen et al. (20041. It is also of interest to compare the peak kinematic growth rates to 
the scaled Lyapunov exponents that are shown in Figure [3jb) . In all cases, it is clear that 
the growth rates are significantly smaller than the corresponding Lyapunov exponents. This 
is unsurprising given that the magnetic Reynolds numbers are comparatively modest in these 
numerical simulations, and we would expect to see higher growth rates at higher values of 
Rm ■ Nevertheless our results would appear to be consistent with the findings of Brandenburg] 



et al. Q1995J) who, for a related dynamo calculation, found that the Lyapunov exponents were 



systematically larger than the observed kinematic growth rates. 

4.2 Magnetic fields in the nonlinear regime 

We now consider the fully nonlinear set of governing equations (Equations |TJ to |6|) including 
the back-reaction of the magnetic field on the velocity field. For each of the four cases under 
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Figure 6: The time evolution of the normalised mean kinetic and mean magnetic energy densities 
for (a) 9 — 3 and (b) 9 — 10. Both quantities are normalised by the mean kinetic energy density 
during the saturated phase (i.e. from t = 200 to t = 500 in the 9 — 3 cases, and from t = 100 to 
t = 250 in the 9 = 10 cases). The magnetic energy has been multiplied by 5. 

consideration, we choose the values of £o corresponding to the largest values of the magnetic 
Reynolds number, in order to maximise the growth rate of the dynamo, thus minimising the 
duration of the kinematic phase. In the following, the time t = corresponds to the insertion 
time of the small magnetic perturbation. The initial ratio between the total magnetic energy in 
the seed field and the total kinetic energy in the flow is roughly 10~ 3 in all cases. This implies 
that the initial field is weak enough to be kinematic, without unnecessarily extending the 
kinematic phase in these nonlinear calculations. Each nonlinear calculation has been evolved 
over a significant fraction of the ohmic decay time (based upon the magnetic diffusivity and 
the horizontal integral scale). 

Figure [6] shows the evolution of the mean kinetic and mean magnetic energy densities for 
the four cases. All of these dynamos are highly time-dependent, but there are clear patterns 
of behaviour. During the early stages of evolution, the magnetic perturbation grows. How- 
ever, the kinematic phase of the dynamo is extremely brief in these calculations, as the seed 
field rapidly becomes dynamically significant. After the short kinematic phase, the dynamo 
undergoes a more extended period of nonlinear growth, eventually settling down to a time- 
dependent saturated state. Conservation of energy implies that the kinetic energy decreases 
during the nonlinear phase of the dynamo, reaching a final level that is clearly lower than the 
kinetic energy of the initial hydrodynamic state. Across all simulations, the mean magnetic 
energy density, (|B| 2 ) /2, saturates at a level that is somewhere between 4 and 9% of the 
mean kinetic energy density, (p|u| 2 ) /2. Although not shown here, a corresponding nonlinear 
calculation for the Rib case (described in the previous section) saturates at a similar level. 
These nonlinear results suggest that the saturation level of the dynamo is comparatively in- 
sensitive to the level of stratification within the domain. There is weak evidence to suggest 
that the rotating cases are saturating at a slightly higher level (particularly for 9 — 3), despite 
the fact that the mid-layer values of the local magnetic Reynolds number are actually smaller 
in the rotating calculations than they are in the corresponding non-rotating cases. If we were 
comparing nonlinear calculations at similar values of the (mid-layer) local magnetic Reynolds 
number, as opposed to similar values of the global Rm, we would expect the rotating cases to 
saturate at a higher level. However, we did not investigate this parameter regime here given 
that this would require higher spatial resolution in the rotating cases. 

Looking again at Figure [6] it is worth noting that the time-dependence of the magnetic 
energy in the non-rotating = 3 case is strongly intermittent. This is clearly seen at t « 
270, where a burst of magnetic energy corresponds to a drop in the kinetic energy. These 
fluctuations in the magnetic energy make it more difficult to determine the level at which the 
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Figure 7: (a) The horizontal magnetic energy, normalised by the equipartition energy B\ Jz) = 
(p|tt| 2 ) z . (b) The equivalent plot for the vertical magnetic energy density, (c) The (unnormalised) 
horizontal magnetic energy density, (d) A plot of B\ q (z) as a function of depth. Note that figures 
(c) and (d) are presented in semi-log scale. 



dynamo saturates in this case. As discussed in the previous subsection, this non-rotating = 3 
case may have a slightly higher critical global magnetic Reynolds number than the other three 
cases. If the dynamo is indeed closer to criticality than the other cases, we would expect it to 
be more sensitive to time-dependent variations in the flow. This would explain the observed 
behaviour. 

A more detailed description of the magnetic field in the saturated phase can be obtained 
by considering the horizontal and vertical components of the magnetic energy density. When 
comparing different cases, it is useful to normalise the magnetic energy densities by the local 
equipartition energy defined by i?e q (z) = (p|it| 2 ) z . The normalised horizontal magnetic en- 
ergy density (B^.) + (By} and the normalised vertical magnetic energy density (-Bf) are 
shown in Figures U[ a) and (b). The unnormalised horizontal magnetic energy density is shown 
in Figure [T^c), whilst the depth dependence of the equipartition energy B% q (z) is shown in Fig- 
ure [7|d). In the middle of the layer, the horizontal magnetic energy density is comparable to 
the vertical magnetic energy density, whereas it clearly dominates close to the boundaries, as 
expected from the chosen boundary conditions. It is clear that the (unnormalised) horizontal 
magnetic energy density is stronger at the lower boundary than it is at the upper boundary, 
although these are more comparable when they are scaled in terms of the equipartition field 
strength. Note that all of these magnetic energy densities are sub-equipartition. Although 
there are some quantitative differences in these curves, the variation with depth of each of 
these quantities is broadly similar in most cases. The main difference to note is that the 
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Figure 8: Grey-scale plots of the horizontal component B x of the magnetic field at two different 
depths, z — 0.2 and z = 0.8, for both rotating and non-rotating cases in the saturated phase. The 
thermal stratification is 9 = 3 in all cases. Contours are evenly spaced between B x = —0.5 and 
B x = 0.5 for Ta = and between B x = —0.3 and B x — 0.3 for Ta = 10 5 . Light and dark tones 
correspond to opposite polarities. 
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Figure 9: The horizontally-averaged ohmic dissipation, (qk (j 2 ) , for (a) 9 = 3 and (b) 8 = 10. 



vertical component of the magnetic energy density tends to be larger in the rotating cases. 

Figure [8] shows some examples of contours of the horizontal component of the magnetic 
field, B x , during the saturated phase. In all cases, 9 = 3, and two different depths are 
considered (z = 0.2 and z = 0.8). The upper row corresponds to the non-rotating case, Rl, 
whereas the lower row corresponds to the rotating case, R2. It is useful to compare this 
figure with the corresponding temperature and enstrophy plots in Figure [3] Near the upper 
boundary (i.e. at z = 0.2), there is a clear correlation between the convective downflows and 
the strongest magnetic fields. The characteristic scale of the magnetic field is smaller in the 
rotating case and one observes that the helical structures of the fluid motion (see Figure [3| 
affect the magnetic field. Close to the lower boundary (i.e. at z = 0.8), both cases are similar 
even if the characteristic scale of the magnetic structures is still smaller in the rotating case. 
Note also that the same colour table is used for both upper and lower layers. The magnetic 
field is therefore stronger and more intermittent near the lower boundary in both cases, as 
already illustrated by Figure [7]^c) . 

The magnetic field geometry plays a crucial role in determining the rate at which magnetic 
energy is dissipated via ohmic diffusion. This is clearly an important process to consider 
when comparing the efficiency of convectively-driven dynamos. In this model, the rate of 
magnetic energy dissipation (per unit volume) due to ohmic diffusion is given by Co K |j| 2 , 
where j = V X B is the current density. This expression for the ohmic dissipation can be 
derived by taking the scalar product of the induction equation Q with B, and then integrating 
over the domain. Note that the rate of dissipation is also proportional to the ohmic heating 
term in Equation Q. The rate of ohmic dissipation in the nonlinear phase, averaged over time 
and horizontal coordinates, is plotted in Figure [jja) for 9 = 3 and Fi gure[9|b) for 9 = 10. The 
dissipation is much larger in the 9 = 10 case, since both the kinetic and the magnetic energy 
densities are greater by roughly one order of magnitude. For each value of 9, the magnitude 
of the dissipation term near the top of the layer is roughly the same, whether or not the layer 
is rotating. On the other hand, rotation clearly reduces the dissipation of magnetic energy 
near the lower boundary. Put another way, the dissipation term is more weakly dependent 
upon depth in the rotating cases than it is in the corresponding non-rotating calculations. 
A similar reduction of the magnetic dissipation by rotation has already been observed in 
rotating homogeneous magnetohydrodynamic turbulence by Favier et al. (20111. It is also 
worth noting that this reduction in magnetic dissipation is not a nonlinear effect: it can 
also be observed in the kinematic phase (although it is more difficult to quantify such a 
reduction since it becomes necessary to normalise the dissipation term by an exponentially 
growing field in order to take an appropriate time-average). The regions of strongest ohmic 
dissipation in the non-rotating cases coincide with the regions of strongest shear. We have 
already seen, in Figures [3ja) and|3|b), that the stretching is more homogeneously distributed 
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Figure 10: Energy spectra for 9 = 3 (the results are similar for 9 = 10). The solid lines correspond 
to the saturated phase whereas the dotted lines correspond to the kinematic phase. The kinetic 
energy spectrum EK(k±_) is averaged over time and vertical coordinate. The magnetic energy 
spectrum is normalised by the time-average of the magnetic energy for the saturated dynamo. 



across the layer in the rotating calculations. As a result, the dynamo-generated magnetic field 
organises itself in such a way that the rate of dissipation in these cases is lower than it is in 
the equivalent non-rotating calculations. So even though the scaled Lyapunov exponents, as 
shown in Figure [3jb), suggest that there is generally less stretching in the rotating cases, this 
reduction in dissipation explains why rotating convection can still act as an efficient dynamo 
(at least at moderate magnetic Reynolds numbers). 

4.3 Saturation mechanisms 

In the previous subsection, we discussed nonlinear results, without really saying anything 
about the saturation mechanism for these dynamo calculations. Due to the complexity of these 
dynamo models, it is difficult to say anything definitive about this. However, some insights 
into the saturation process can be gained by comparing certain aspects of the kinematic and 
nonlinear phases of the dynamo. 

Firstly, we consider some of the global properties of the dynamo. We define the magnetic 
energy spectrum in the following way: 

= ^^^B^ky^) • B*(k x ,k y ,z) ( 18 ) 

2 kj_ 

where k\ — k x + k y is the horizontal wave number, B(k x , k y , z) is the two-dimensional Fourier 
transform of B(x,y,z) and the star denotes the complex conjugation. Similarly, the kinetic 
energy spectrum is defined as follows: 

E K (k±) = ~ y^y^ u(k x ,k y , z)-pu*(k x , k y ,z) + u*(k x ,k y , z)-pu(k x , k y , z) . (19) 

In Figure |10| we show the energy spectra corresponding to the 8 = 3 calculations (Rl and 
7?2), in both the kinematic and nonlinear regimes. The kinetic energy spectrum is time- 
averaged, whereas the magnetic energy spectrum is normalised by the time-average of the 
magnetic energy for the saturated dynamo before also being time-averaged. The results are 
qualitatively similar for the highly-stratified case. Comparing the kinetic energy spectra, we see 
that there is less kinetic energy at large scales (i.e. for 1 < k± < 3) in the rotating case. This 
is connected to the rotationally-induced reduction in the horizontal scale of motion. Whether 
or not rotation is present, the magnetic energy spectra always peak at small scales, which is 
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Figure 11: Horizontally-averaged maximum Lyapunov exponents versus depth during the nonlinear 
phase (solid lines). The dotted lines correspond to the previous kinematic results (see figure [31(a)). 



consistent with the fact that we do not observe a large-scale dynamo in these simulations. 
Comparing the kinematic and nonlinear phases, we see that the kinetic energy spectra are 
only weakly perturbed by the magnetic fields. So there is no evidence here to suggest that 
saturation is accompanied by a significant variation in the kinetic energy spectrum. More 
interestingly, there is a small (but perhaps significant) alteration of the magnetic energy spectra 
as the dynamo saturates. There is a clear reduction of magnetic energy at small scales and 
a corresponding increase at large scales. This is observed in both the rotating and the non- 
rotating cases. Although not shown here, the same trend is observed in the — 10 calculations. 
We note that this is consistent with the increase of the Taylor microscale of the magnetic field 



reported by Brandenburg et al. (1996) 



When considering potential saturation mechanisms, one of the most interesting things to 
consider is whether or not the stretching properties of the flow are modified by the magnetic 
fields. Local variations in the stretching would probably not lead to significant variations in the 
kinetic energy spectrum. However, we would expect to see changes in the Lyapunov exponents 
if the stretching properties of the flow are modified in the nonlinear dynamo regime. Using the 
flows from our nonlinear calculations, we evaluate the Lyapunov exponents using exactly the 
same procedure as for the kinematic phase (as described in Section |3.2[ ). The particle pairs are 
followed from t = 300 to t = 500 in the — 3 cases, and from t = 150 to t — 250 in the = 10 
cases. The (horizontally-averaged) maximum Lyapunov exponents for the four cases are shown 
in Figure [TT] by the solid lines. The dotted lines show the corresponding Lyapunov exponents 
from the kinematic calculations. In the non-rotating cases, regardless of the value of 0, the 
Lyapunov exponents are slightly lower in the nonlinear phase everywhere except near the upper 
and lower boundaries. In the rotating cases, the Lyapunov exponent is slightly reduced almost 
everywhere, even close to the lower boundary. So there is some indication of a suppression of 
stretching due the presence of magnetic fields. Furthermore, given the intermittent nature of 
the magnetic field distribution, we might expect the local reduction in stretching (in regions 
of strong magnetic fields) to be greater than that suggested by these horizontally-averaged 
quantities. However, this does not mean that the saturation of the dynamo can be explained 



simply by a reduction in stretching. Cattaneo & Tobias ( 2009 1 have shown that a dynamo- 



saturated velocity field in Boussinesq convection can still act as a kinematic dynamo (provided 
that the new seed field is not aligned with the original magnetic field). Although this has not 
been tested here, the depth-dependence of the Lyapunov exponents in the nonlinear regime is 
qualitatively similar to that of the original hydrodynamic flows, so we would speculate that 
these dynamos would exhibit similar behaviour. Certainly we do not see the drastic reduction 



in the Lyapunov exponents that Cattaneo et al. ( 1996 \ observed in their model of nonlinear 



dynamo action in a much simpler one-scale flow. The saturation process in convectively-driven 
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Figure 12: Alignment properties for the simulation R2 (the results are qualitatively similar for our 
other simulations), (a) Probability density function of the cosine of the angle between u and B. 
We separate contributions from points where \B\ > 3£> rms (solid lines) and \B\ < 3B lms (dotted 
lines), (b) Same as (a) but for the angle between B and e 3 , where e 3 is the eigenvector associated 
to the largest eigenvalue of the rate of strain tensor. 



dynamos appears to be more subtle than this. 

We also consider the possibility that the correlation between u and B plays some role in 
the saturation process. The alignment between these two vectors, and the influence that this 
has upon some of the key nonlinearities in the dynamo system, has been extensively studied 
in recent years (see, for example, Servidio et al. 20081. Figure 12 a) shows the probability 
density function of the cosine of the angle between u and B: 



cos(u, B) — 



\u\\B\ 



(20) 



The pdfs have been obtained by averaging over space (between z = 0.2 and z — 0.8 only, so 
as to focus upon the interior of the domain) and time (using around 50 snapshots in each 
case). Note that it is not necessary to carry out any normalisation to compare pdfs in the 
kinematic and nonlinear regimes. Given that we would expect any preferential alignment to be 
rather localised, we separate regions of strong magnetic fields from weaker field regions. This 
is achieved by plotting two pdfs, one for those mesh points for which \B\ > 3-B r ms, where firms 
is the rms magnetic field, the other for those mesh points that fall below this threshold field 
strength. Despite the fact that this flow is compressible, inhomogeneous and anisotropic, we 
see remarkable alignment between u and B. Perhaps unsurprisingly, the alignment is always 
slightly more pronounced for strong magnetic fields than for the weaker fields. However, the 
strong alignment between the two fields is observed in both the kinematic and the nonlinear 
phases, and is therefore a property of the induction equation. Hence the saturation of these 
convectively-driven dynamos cannot be explained by a modified alignment between u and B. 



Following Brandenburg et al. ( 1996[), we also consider the alignment between B and the 



eigenvectors of the rate of strain tensor S. The symmetric matrix Sij = | (uij + Uj ti ) is com- 
puted at each mesh point and the three corresponding eigenvalues are ordered from the smallest 
to the largest. The smallest eigenvalue is always negative and its eigenvector corresponds to 
the direction of compression. The magnetic field is found to be mostly perpendicular to this 
direction (as already reported by Brandenburg et al. 19961, in both the kinematic and the 



nonlinear phases. The largest eigenvalue is always positive and its eigenvector, denoted here 
by e3, corresponds to the direction of maximum stretching. The probability density function 
of cos(B,e3), is shown in Figure 12 b) . Note that we again separate the contributions from 
the weak field and strong field regions. During the kinematic phase, the strongest magnetic 
fields are mostly aligned with the direction of maximum stretching. Note that there is also 
some indication of preferential alignment for the weaker field regions, but this is much less 
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pronounced. In the saturated phase, the pdf corresponding to the weak field regions shows 
a very modest reduction in alignment. However, there is a dramatic reduction in the align- 
ment between B and e$ in the strong field regions. Note that the magnetic fields for which 
\B\ > 3-B rms only represent 2% of the total number of points, so that the loss of alignment 
between the direction of maximum stretching and the magnetic field is only observed very lo- 
cally. This modification of the alignment between B and is observed in all our simulations, 
with or without rotation, and for both thermal stratifications. This indicates that this process 
is rather robust. Given that this modified alignment will reduce the efficiency of the dynamo 
in regions of strong magnetic fields, we conclude that this effect plays a role in the saturation 
of convectively-driven dynamos. 



5 Conclusions 

In this paper, we have investigated small-scale dynamo action in rotating compressible convec- 
tion. Regardless of the level of stratification within the domain, both rotating and non-rotating 
convective flows can sustain a small-scale dynamo if the magnetic diffusivity is small enough. 
Using a mid-layer value for the magnetic Reynolds number (based upon the integral scale of 
the turbulence rather than the layer depth), rotation seems to reduce the value of the magnetic 
Reynolds number above which dynamo action is observed. Increasing the thermal stratifica- 
tion also reduces the critical value of the local magnetic Reynolds number in the non-rotating 
case. At high values of the magnetic Reynolds number, the growth rate of the magnetic energy 
of the dynamo appears to have a logarithmic dependence upon Rm in the weakly-stratified, 

1/2 

non-rotating simulation. It is more difficult to fit a scaling law in the other cases, but an R AI 
scaling may be more appropriate here. At a given value of the mid-layer (local) magnetic 
Reynolds number, the kinematic growth rates are always larger in rotating convection than 
they are in the corresponding non-rotating cases. This dependence upon rotation cannot be 
attributed solely to the fact that the local Reynolds number is smaller in the rotating cases. 
So we conclude that the Coriolis force plays a key role in determining the kinematic dynamo 
properties of rotating convection. At the highest computationally accessible values of the mag- 
netic Reynolds number, we find similar levels of saturation in all of our nonlinear calculations 
(with the magnetic energy saturating at about 4 — 9% of the final kinetic energy). At first 
sight, this result is slightly surprising given that the Lyapunov exponents suggest that there 
is less stretching in the rotating cases (particularly near the lower boundary). However, this 
is compensated by the fact that magnetic dissipation seems to be much less efficient in the 
rotating calculations. It is difficult to say anything definitive about the saturation mecha- 
nism for these convectively-driven dynamos. However, a comparison between the kinematic 
dynamo regimes and the nonlinear saturated states show a slight reduction in the Lyapunov 
exponents in the nonlinear regime (due to a local reduction of the stretching properties of the 
flow). Furthermore, there is some evidence to suggest that a reduction in alignment between 
the strongest magnetic fields and the direction of maximum stretching also plays a role in the 
saturation process. 

The next natural step in this study is to explore the parametric dependence of this system 
in more detail. Of course, given the complexity of this numerical model, it is not feasible 
to conduct a complete survey of parameter space, although any further calculations that can 
be carried out will obviously enhance our understanding of this problem. We also intend to 
investigate the ability of rotating compressible convection to produce a magnetic field on much 
larger scales than the scales that are associated with the convective motions. The question of 
how hydromagnetic dynamos are able to generate large-scale magnetic fields is undoubtedly 
one of the most challenging issues in modern dynamo theory. In the standard formulation of 



mean- field dynamo theory (Moffatt 19781, this regenerative process relies upon the presence 
of helical motions, which are invariably produced when a flow is influenced by rotation. Hence, 
rotating convection should be able to drive a large-scale dynamo. However, despite predictions 



from mean-field dynamo theory, Cattaneo & Hughes ( 2006 ) failed to find evidence for a large 



scale dynamo in their Boussinesq model. Interestingly, Kapyla, et al. (2009) did find large-scale 
magnetic fields in their compressible model. They argue that the absence of large-scale fields 
in the Boussinesq model of Cattaneo and Hughes can be attributed to a rotation rate that is 
too slow. This is a possibility, however the effects of compressibility may also be playing a 
role. The flow is (locally) strongly helical in the Boussinesq calculations of|Cattaneo fc Hughes| 



(2006), but the mid-layer symmetry of their setup implies that the mean helicity distribution 



is antisymmetric about the mid-plane. It may be important to break this symmetry in order 
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to generate large-scale fields. In work that is currently in progress, we are investigating this 
issue by carrying out simulations of rapidly-rotating convection in a wide compressible layer. 
The size of the layer is important, because it is necessary to have a clear separation in scales 
between the small-scale fields and any large scale magnetic fields that may be generated. The 
computational domain that was considered in the present study was too small to allow for this 
separation in scales, which may explain why no large scale fields were observed here (despite 
the fact that the helicity of the flow is asymmetric about the mid-plane in these stratified 
rotating calculations). Of course, other physical ingredients could also be included into this 
system once the basic ingredients of rotation and stratification are properly understood. For 
example, the role of the penetrative layer in the model of |Kapyla et al.| (|2009) is also unclear. 
This may be promoting the formation of large-scale fields in some way. It may also be of 



interest to consider the effect of including a shear flow (see, for example, Hughes & Proctor 
2009) in this compressible model. 
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